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Summary 

The accuracy is computed of a computational procedure for rapid estimation of frequency response 
functions from stereoscopic dynamic data. It is shown that reversal of the order of the operations of coor- 
dinate transformation and Fourier transformation, which provides a significant increase in computational 
speed, introduces error. A portion of the error can be made arbitrarily small by prescaling frequency 
data prior to coordinate transformation. The remainder of the error, which cannot be eliminated, results 
from components of motion normal to the camera focal planes. Spectral analysis of the dominant error 
term shows the presence of harmonic and cross-modulated spectral peaks caused by the nonlinearity of 
the perspective projection from spatial coordinates onto focal plane coordinates. Only when the camera 
geometry is such that all motion is parallel to both camera planes is the spatial transformation linear and 
no error introduced. 

The full 3x3 frequency response matrix of a vibrating structure can be estimated from a single stereo- 
scopic dynamic data set only if the driving function consists of three orthogonal mutually independent 
stochastic processes. Otherwise three independent sets of separate measurements are necessary. A least- 
squares procedure for estimation of the full frequency response matrix, column by column, for the latter 
case is developed. 

Introduction 

Close-range photogrammetric techniques have application to remote measurement of displacement 
in aerospace structural testing. Through the use of these techniques, the motion of multiple targets 
on a vibrating surface can be recorded by means of a stereoscopic digital camera system as a pair 
of two-dimensional (or equivalently, a four-dimensional) sampled-data time series. Originally the four- 
dimensional coordinates of each stereoscopic data point, recorded for each target at each sampling time, 
were transformed into three-dimensional spatial coordinates by means of the direct linear transformation 
(DLT) algorithm (ref. 1). Cross power spectral matrices were then estimated in spatial coordinates 
through the use of averaged discrete Fourier transformation (DFT). Kroen and Tripp (ref. 2) showed 
that the operations of coordinate transformation (DLT) and frequency transformation (DFT) may be 
commuted subject to certain restrictions, i.e., limitation of displacement normal to camera focal planes. 
The computation time per point required for Fourier transformation and averaging is negligible compared 
with that for coordinate transformation, which is required only for each averaged data point in the 
DFT-DLT transformational sequence. Thus, reversal of the order of operations produces a significant 
computational time savings (typically a factor of 125, ref. 2) at the cost, however, of error due to the 
nonlinearity of the direct linear transformation whenever components of motion exist normal to either 
camera focal plane. 

This paper derives an analytical upper bound on the magnitude of the error incurred by the DFT-DLT 
operation sequence as a function of test geometry and perturbation amplitude. A procedure using the 
DFT-DLT transformational sequence is also developed for least-squares estimation of the general 3x3 
frequency response matrix. 

Symbols 

A projection matrix defined in equation (20) 

Aj. j th row of matrix A 

a damping rate; power series coefficient; general scalar element 

(iij element of matrix A 

a\ defined in equation (111) 

B general matrix 



bi 

bnm{t) 
b A 

C 

C(6w) 

C[F(<5w)] 

c 

j 

' mk- ' ink 

D 
Da 
D z (t) 
d 

dj(t) 

d X 

E 

E(w) 
EkM 
e r(f) 

'V(0 

e A '(0 

/'[/(')] 

f(t) 

&(t) 

H(u) 

h 

ba(tj) 

hz 

I 

i 

j\ m. 71 , r, q 

K 

M 

M 1 ’ 2 


ratio defined in equation (107) 
power series coefficient 
defined in equation (112) 

maximum perturbation magnitude defined in equation (74) 

perturbation matrix defined in equation (24) 

matrix defined in equation (50) 

camera focal length 

power series coefficients 

mapping defined in equation (28) 

mapping defined in equation (21) 

diagonal matrix defined in equation (70) 

power series coefficient 

power series defined in equation (118) 

defined in equation (113) 

expected value 

error defined in equation (35) 
error defined in equation (03) 
error defined in equation (05) 
defined in equation ( 1 10) 
error defined in equation (00) 

Fourier transformation of function f(t) 

perturbation function 

driving function 

frequency response matrix 

vector in OBJ coordinates 

zth diagonal element of H(cj) 

defined in equation (157) 

j th column of H(u;) 

defined in equation (58) 

identity matrix 

V^-T: integer 

integers 

scale factor 

rotation matrix defined in equation (1) 
rows 1 and 2 of mat rix M 
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M 3 

M ff 

N 

Pa 

Pn 

R uuW 

Rwu( r ) 

Suu( w )> Suu(w) 

SwuM) Swu M 

s ggi UJ ) 

Suu.j (^) 

T 

t 

U* 

6U{t) 

fiU(w) 

u,v 

<5u(£) 

W/ 

W 

w 7 

«W(t) 

«W(w) 

w 

<5w(£) 

x c . y c , z c 

Xp,Yp, Zp 

x.y.z 

Z 

«z 

«Z f [F(«w)] 
a, 0 

Y k 


row 3 of matrix M 

composite matrix defined in equation (13) 
number of averaged records 

generalized inverse of matrix A defined in equation (23) 

constant matrix 

autocovariance matrix 

cross covariance matrix 

autospectral matrices 

cross spectral matrices 

autopower spectrum of g(t) 

j th column of S U u(w) 

record length, seconds 

time, seconds 

constant vector 

input driving vector in object space 
= F[SU(t)} 
general vectors 

input driving vector in DFP space 

zth element of v 

point in object space 

coordinates of point W 

focal plane image of W 

normalized magnitude of 6W(t) 

time- varying perturbation in object space 

= F[6W{t)] 

point in DFP space 

time-varying perturbation in DFP space 

object space coordinates 

object space coordinates of perspective center 

focal plane coordinates 

diagonal matrix defined in equation (19a) 

diagonal perturbation matrix defined in equation (19b) 

diagonal matrix defined in equation (40) 

constants 

matrix defined in equations (162), (165), and (167) 
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Subscripts: 

k 

P 

1 


2 


Superscripts: 

T 


* 


least value of ^(t) 
element of A^ 

angle between 6W(£) and ^-axis 

perturbation symbol 

RMS error defined in equation (173) 

maximum relative perturbation defined in equation (73) 

relative RMS error defined in equation (82) 

RMS error defined in equation (171) 

maximum negative relative perturbation defined in equation (78) 
vector defined in equation (11) 
orientation angles (fig. 3) 

= Z _1 Ar fc 

eigenvalue of (A T A)~ l 
jth column of A*, 
constant defined in equation (53) 
perspective center coordinates 
transformed autospectral matrix 
transformed cross spectral matrix 

condition number of matrix A defined in equation (69) 
elements of £ uu (u;) and £ wu (u;) 
jth column of £ wu (u;) 
time 

PC coordinates of point W 
focal plane image of T 
frequency, rad/sec 
damped frequency, rad/sec 

k — 1 for camera 1; k = 2 for camera 2 
perspective center 
camera 1 
camera 2 
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matrix transpost' 
complex conjugate 


Abbreviations: 


DFP 

DFT 

DLT 

FP 

max 

OBJ 

PC 

RMS 


dual focal plane 

discrete Fourier transformation 

direct linear transformation 

focal plane 

maximum 

object 

perspective center 
root mean square 


Mathematical notation: 

= is defined to be 

|| || vector or matrix norm 

A caret (~) over a symbol denotes an estimated value. 


Development of Projection Equations 

The three-dimensional position of a vibrating target point W is remotely sensed as a function of 
time by means of a stereoscopic camera system which employs two cameras focused on W , as shown 
in figure 1. The image of W , denoted by W ! , appearing in each camera focal plane is modeled as a 



W 

object 

Figure 1. Stereoscopic camera geometry showing images of object W on camera focal planes. 


perspective projection. That is, W projects onto the point of intersection of the camera focal plane with 
the line passing from W through the persective center, W p (the lens center, ref. 3). For the perspective 
projection geometry, as shown in figure 2, define three sets of spatial coordinates: object (OBJ) coordinates 
X c ,Yc,Z c ; perspective center (PC) coordinates whose origin is the perspective center; and focal 

plane (FP) coordinates x,y,z whose origin lies on the focal plane. The PC and FP coordinate systems 
are parallel and displaced from one another by focal length c in the ^-direction, but they may be rotated 
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OBJ coordinates c PC coordinates FP coordinates 


Figure 2. Perspective projection coordinate systems showing projection of object W through perspective center onto 
image point W ' in focal plane. 



Figure 3. Angles of exterior orientation showing rotation from coordinates (X C ,Y C , Z c ) into coordinates (x^, y/c , Zk)- 

in space relative to the OBJ coordinate system. The rotational transformation from OBJ orientation to 
PC-FP orientation is computed via rotation matrix M (ref. 3) where 

cos (/> cos 0 sin t/> sin (j) cos 0 4- cos ip sin 0 sin sin 0 — cos sin (f) cos 0 

M= — cos0sintf cos t/»cos 0 — sin 0 sin <j) sin 0 sin ifc cos 0 + cos ij) sin (j) sin 0 (1) 

si \\(f) — sin 0 cos 0 cos cos (j) 


See figure 3 for angle and coordinate definitions. 



Refer to figure 2 for the following steps. Let the object (OBJ) coordinates of object point W be denoted 


by 

W = [X C Y c Z c f OBJ 

and let the perspective center (PC) coordinates of W be denoted by 

T = K V dpc 

Let W p denote the coordinates of the perspective center in OBJ coordinates as 

W p = [X p Y p Z p } 0 B j 

The transformation from OBJ coordinates to PC coordinates is then 


(2) 

(3) 

(4) 
(3) 

OBJ 

Image point Y' — [£' rj f $ f ] T , the point of intersection of the focal plane with the line passing from 
T through the perspective center (the origin in PC coordinates), is obtained by scaling T so that its $ 
component equals — c. Thus, 


T = < 

m 

V 

> = M < 

\x c -x p ) 
Y c ~ Yp 


LJ 

1 PC ■ 

[ Z c — Zp J 


T' = 



( 6 ) 


OBJ 


is the perspective projection point of W in PC coordinates. It is transformed into focal plane (FP) 
coordinates (whose axes are parallel to the PC coordinate axes) by translation to the focal plane. Let 
(zp, 2 / p ,c)fp denote the FP coordinates of the perspective center and let (x,2/,0)pp denote the FP 
coordinates of the perspective projection image point. Then the PC coordinates of the perspective 
projection image point may be expressed as 




PC 



(7) 


Thus the transformation from OBJ coordinates to FP coordinates, obtained from equations (6) and (7), 
is 


( x - Xp ' 

c 

[X c -x p ) 

\v~Vp 

= --M 

Y c -Yp 

\ -c J 

I FP ' 1 

{ Z c Zp J 


( 8 ) 


OBJ 


By means of the two cameras, a subset of the three-dimensional object space is mapped into the two 
focal plane coordinate systems, which can be viewed as a four-dimensional, dual- focal-plane (DFP) space. 
Let subscripts 1 and 2 in the sequel denote cameras 1 and 2. Let M 1,2 denote the 2x3 matrix consisting 
of rows 1 and 2 of M, and let M 3 denote the third row of M. Also let 


and 


w = [xi m 

X2 V2\ T 

(9) 

wi = [x\ 

y i] T 

(10a) 

W2 = [X2 

V2] T 

(10b) 

(m = ki 

f2] T 

(11) 
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Apply equation (5) to cameras 1 and 2 and combine to obtain an expression for $\j as 


Cm = 


(?) 

= M f < 

[£}- 

Mf 

1 

! 0 

lC2 J 


[Zcf 

0 

! Ml 


(Xpi 

Y Pi 

Zp i 
< — 


X 


P2 

Y P2 

Zpi 


( 12 ) 


where 


M c = 


Mf 


(13) 


To obtain the elements of w, similarly apply equation (8) to cameras 1 and 2. Extracting rows 1 and 2 
of equation (8) for each camera and combining gives 


/ *1 - x Pi \ 

l 2 /i- y Pl J 


X\ - x Pl 

y\- ypi 

X2 - X P2 

2/2 - y P 2 


Cl 

C2 


> = 


-ci iMj ; 2 

L-c 2 mJ’ 2 


X,. 


+ 


qM}' 2 

1 

! o 

0 

c 2 M l 2 


X 


Pi 

K P1 

? p i 


X, 


Y, 


P 2 


v z 


P 2 
V 2 


(14) 


Let />W = [<*>A' ( ; 6Y C bZ c ] T lie a perturbation in W. FVom equation (12) the corresponding perturbation 
^C M 18 

' t)X r 


/’Ciw 


i / *^1 1 

/> C2 j 


= M f { 6Y C }■ = M f 6 W 
6Z C 


(15) 


Now substitute the elements of w + bw. Cm + <5Cm- and W + b W into equation (14) for the elements of 
w,c. and W, and subtract equation (14) from the result. This yields 


/>Cl ■ 

' Xi - x Pl ' 

?/i - y P i J 

' + (Cl+«Cl) < 

bxi ) 

byi J 

( i) C2 j 

f •*•2 - X P2 1 

l 2/2 - 2/p2 J 

j + (C2 + 6'C 2 ) | 

f fix* ] 
[ t>V2 j 


> = 


-cM 


1,2 


-cM 


1,2 
2 J 



(16) 


Equations (15) and (16) are combined to obtain 


1.2 


■qMj 
-c 2 M’’ 2 j 


X'2 - X P2 

2/2 - y V2 J 


m]j 


Vi - Vpi 

0 
0 

I - '' r & _i 

or in matrix notation, by lett ing Cw = [Am i by i Am 2 />y 2 ] r 

A tfW = (Z + tfZ)«w 



(ft+*i){J£; 


(17) 


(18) 
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Figure 4. Convergent projection geometry showing projection of object perturbation 6 W onto image perturbations 6 v/\ 
and 6 w 2 in focal planes 1 and 2. 


where 


and 



'SI 

0 

0 

0 - 



0 

Si 

0 

0 




0 

0 

S2 

0 



.0 

0 

0 

C2 - 



r<Hi 

0 

0 

0 ■ 



0 

<5si 

0 

0 



0 

0 

<$S2 

0 



0 

0 

0 

<5?2- 


A = 



2-1 *^Pl 0 

y\ ~ y P i 0 




— C 2 M 2 2 


0 2/2 Xp2 

0 2/2 - y P2 . 


_ M 2. 


(19a) 


(19b) 


( 20 ) 


An illustration of perturbations 8 W and <5w appears in figure 4. 

The projection equation for <Sw as a function of <5 W can be obtained from equation (18) as a mapping 
D a from OBJ coordinates into DFP coordinates as 


6w = (Z + 6Z)~ 1 A6W = D a {8W) (21) 

Matrix (Z + 8Z)~ l exists whenever no perturbation in the ^-direction intersects the (£, 77 ) plane (which 
intersects the perspective center and is parallel to the focal plane, as shown in fig. 2). Equation (21) 
defines 8 w as a function of 8 W, since the elements of 8 Z are obtained from equation (15), which is 
dependent only on <5W. If matrix A has full rank (rank 2), then from equation (18) one can obtain 

8W = P A (Z + <5Z) tfw (22) 

where 

P A — (A^Aj 1 A t (23) 


9 








Matrix Pa, a generalized inverse (ref. 4), is the two-sided inverse of A over the range (column space) 

of A. 

It is now shown that 6 Z is dependent only on 8w in equation (22) and that equation (22) therefore 
defines a mapping D from DFP coordinates into OBJ coordinates which is the inverse of mapping D 4. 
Define a matrix containing the elements of 6 w (with the help of eq. (9)) as 

~6x\ 0 ' 

C(Sw)= s l 2 (24) 

- 0 8t)2 . 

which yields the identity 

6 Z <5w = C(<5w) (25) 

From equations (15). (22), and (25) one obtains an expression for 

S( M = M f P A Z bw + M f P A C(*w) S( M (26) 

which is solved for 8 (m to obtain 

S( M = [I - M f P A C(^w)] -1 M ? P a Z tfw (27) 

Thus, whenever [I - M^Pa C(<$w)] is nonsingular, equations (22), (25), and (27) can be combined to 
obtain 

«W = P a (z + C(M [I-M c P a C(J»w)] _1 M f P A z} <5 w = D{8 w) (28) 

From equations (21), (22), and (28) it. is seen that 

8W — D[D A {t> W)] (29) 

and I) is a left inverse of 1)^. In the applications for which this study is intended, every value of 8w will 
be contained in the range of mapping D A \ that is, will be the image under D A of some perturbation 
dW in OBJ coordinates. Thus, if 6w — ^^(^W), then 

D a [D( 8 »w)] = D A |D[D^(«W)]} = D A (8 W) = tfw (30) 

Honco. I) is 1 ho right inverse of D A over the range of D A . Equation (28) shows that, in general, D is 
nonlinear in iw. 

Error Caused by Commutation of Coordinate Transformation and Fourier 
Transformation 

Let ^W(t) and <Sw(/,) denote time-varying perturbations of <^W and t>w. (The (l) notation will often 
be suppressed in the sequel except where needed for emphasis or clarity.) Fourier transformation of vector 
functions 8W(t) and 8w(t) is obtained by componentwise scalar Fourier transformation as 

( 6X e (u) ) 

F[bW(t)} = l 6Y e (ui) } (31) 

{ bZ r (u) ) 


and 


F[(')w(0] 


^xi(u;) > 

by iM 
bX2{u) 
^2H > 


(32) 
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which are complex-valued vector functions of u. Mapping D is extended to complex-valued vectors by 
the relation 

D{u + iv) = D(u) + i D{v) (33) 

If D were linear, then (since Fourier transformation is a linear operation) composition of mapping under 
D and Fourier transformation would be commutative. However, because of the nonlinearity of Z), it 
follows that 

F [6W{t)] = F {D [tfw(0]} + D {F [Sw{t)]} (34) 

If the perturbation magnitudes are small, the error 

E{uj) = D{Fl6v,(t)}}-F[6W{t)} (35) 

may be acceptably small. An expression for E(u;) is now obtained. 

It is easily shown that for nonsingular matrices Z and Z + SZ 

(Z + <5Z) -1 = Z -1 |l — <5Z(Z + <5Z) _1 ] (36) 

Substitute equation (36) into equation (21) to obtain 

tfw(i) = Z~ l {A 6Vf(t)-6Z(t)[Z + 6Z(t)]~ l A «W(t)} (37) 

The Fourier transform of equation (37) is then 

F(<5w) = Z" 1 {AF(<5W) - F{8Z { Z + 6Z)~ l A 6W]} (38) 

To obtain D[F(6 w)], substitute F(6w) into equation (28) for <5w to yield 

D[F{6vr)] = P A {Z + 6Z f [F{6w)}Z}F{6w) (39) 

where 

<5Z f [F(<5w)] = C[F(6w)] {I - M ? P a C[F(<5w)]} _1 M ? P a (40) 

and C[F(<5w)] is obtained by substitution of F(8 w) for 8 w into equation (24). Combine equations (38) 
and (39) to obtain 

F>[F(«5w)| = F(<5W) - P A F [$Z(Z + 6Z)- 1 A6\v] 

+ P A <5Z F [F(6w)]{AF(6W)-F[<5Z(Z + 6Z) 1 A <5w] J (41) 

The expression for E(cj) is obtained from equations (35) and (41) as 

E(w) = — P a {I + <5Z f [F(<5w)]}F [tfZ(Z + «Z) _1 A$w] + P A <5 Z f [F(<Sw)] A F(fiW) (42) 


Determination of Error Bounds 

Let || v || denote the Euclidean norm of an n-element vector v, where 




1/2 


(43) 
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and let ||A|| denote the norm of a rectangular matrix A, where 


l|A|| 


A 

= max 
v^O 


l|Av|| 

l|v|| 


( 44 ) 


It can be shown that the value of ||A|| is equal to the largest singular value of A (ref. 5). The value ||A *||, 
when A -1 exists, equals the inverse of the smallest singular value of A. It follows from equation (44) 
that 


II Av|| < 1| A || • ||v 


(45) 


and that there always exists a v for which equality holds in equation (45). 
I hat 


IABH < ||A| 


IB I 


It can also be shown (ref. 5) 

(46) 


where A and B are conformable matrices. The triangle inequality for vector norms given by 


||u + v|| < ||u|| + ||v|| (47) 

is true also for matrix norms (ref. 5) as 

II A + B|| < || A || + || B || (48) 


In the sequel, it is assumed, realistically, that the norms ||tfW(£)|| and ||F(/>W)|| are bounded, and that 
the integral \\F(b W)|| 2 dw exists. 

Take the norm of equation (42) and apply inequalities (45) to (47) to obtain 


l|EM||<||P A |M|I + ^Z / ,[/^w) 


F 


6Z(Z + tiZ)- l AfiW 


+ l|P A || • II A|| • \\6Z f [f(iw)l II • ||F(/>W)|| 


(49) 


Inspection of equation (49) shows that when ||/>Z/.'[F(i'w)]|| is sufficiently small, E(u) is dominated by 


Pa F 


<SZ(Z + <SZ)- 1 A<SW 


It is now shown through the use of equation (40) that 

be made arbitrarily small by appropriate scaling of F(6w). 

From equation (24) it is seen that 


>Z /;’[/'■’(<'> w)] || can 


C[F(6w)} 


'Axi(u;) 0 

byiW) 0 

0 fix 2 (cj) 

0 ^?/2( w ) - 


(50) 


It is readily shown that the norm of C[F(/>w)] is bounded by 

||C[F(/>w)]|| = inax[||<Swi(w)||,||<5w2M||] < ||F(«!>w)|| (51) 

It can be shown (ref. 4) that for any square matrix A, if || A|| < 1. then (I — A) -1 exists, and its norm 
is bounded by 

ll(I-Ar ' lls rqjAii (52) 

Let 

" = l|M s -|| • ||P A || ■ ||F(tfw)ll (5») 

It follows from equations (■!(>), (51), and (53) that 

||M s -P A C[F(Fw)]|l < v (54) 
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By the result cited above (inequality (52)), if u < 1, then {I-M^PaC^w)]} 1 exists, and thus 
<5Z^[.F(<5w)] (eq. (40)) exists. It follows from equations (52) and (54) that 


{I — M ? P A C[F(6w)]} 


-1 


< 


3W 


l-||M ? P A C[F(«5w)]|| - 1 - u 

Take the norm of equation (40) and use inequalities (46), (51), and (55) to give for u < 1 

||«Z f [F(tf«;)]|| < ||M f || • ||P A || • II F{6y 

Inequality (56) may be rewritten as 

||<SZ f [F(<5w)]|| < h z 

where 


(55) 


I max 


/ lz = ||M ? ||-||P A ||( r i-)||F(^ 

Let F{6 w) be scaled by factor so that equations (53) and (57) become, respectively, 


^=-||M f ||-||P A || • ||F(£w)||max 


and 


SZ F 


— F(6w) 


< 


K 


(56) 

(57) 

(58) 

(59) 

(60) 


Thus K can be chosen sufficiently large to ensure that v < 1 is satisfied and to make 6Z p \jtFV w) 
arbitrarily small, as was to be shown. 

The effect of scaling F(6w) in equation (49) can now be determined. Replacement of Z?[F(tfw)] in 
equation (39) by KD ^i r ’(<5w)j yields 


KD\-F{6m,) 


{ 


= P A lZ + 6Z F 


— F(<5w) 


i) F(8w) 


Define scaled error E a (cj), corresponding to E(w) in equation (35), as 


K 


F(6vf) 


— F(6W) 


E K {u) = KD 

The expression for E^(w) analogous to equation (42) is written as 

Ejc(w) = -P A |/ + 5Z f ~F{6w) ) P [i5Z(Z + 6Z) _1 Ai5wj 


(61) 


(62) 


+ P A 


—F{6w) 


AF{6W) 


(63) 


Through the use of equations (60) and (63) it follows that the norm ||E^{cj)|| is bounded by 

UEtfMII < (l + ^*) IIPaII • {<5Z(t)[Z + «5Z(t)]- 1 A5W(t)} || 


+ ^*||PaIM|A|M|p[«w(<)]|| 


(64) 
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The first Fourier transform term in equation (64) cannot be analytically evaluated even for simple 
time functions because of the inverse factor [Z + 8Z[t)]~ l . This fact prevents determination of a useful 
bound. However, the integral of ||E#(u;)|| 2 over cj can be bounded analytically, as will now be shown. 
Let 

e F (t) = \\P A \\6Z(t)[Z + 6Z(t))- l A6W(t) (65) 

and 

e*(0 = (l + e F (t) + ^^||P A || • ||A||*W(i) (66) 

From equations (64), (65), and (66), it is seen that 

IIEa'MII < || J F[e /f «)]|| (67) 

Take the norm of equation (65) and apply inequalities (45) and (46) to give 

|MOII<a|lD,(t)|| -||tfW(0|| (68) 

where 

a = ||P A || ■ || A|| (69) 

and 

D z (t) = 6Z(t)[Z + 6Z(t)}~ 1 (70) 

Each nonzero element of D z (t), which is diagonal, is of the form 

(0 A/c 

1 + >k] 

where k — 1 or 2. Factor 8^(1) is obtained from equation (15) as 

*ftfc(0 = M%6W(t) - c(»7*(/.)||tfW(f)|| (71) 

where ')k(t) is the angle between 8W (/) and the ^-Jixis. If % is the smallest angle attained for any t 
between <SW(£) and the ^-axis, then it follows from equation (71) that 


«ft(0| 

< 

cos 



C/c 


>W(* 


and that 

where 

Define 


A 

— max 

<5ft fc(0" 

= C max 

cos ^ 

kd 

C/c 

k 

$ k 


C = max||tfW(0|| 


^f(0 = ^PW(0H 


Combine equations (72) through (75) to obtain 

I %(<) 


St k 


< ( m W F (t) 


Then it follows from equations (70) and (76) that 


1^2 (Oil < ~W F (t) 

1 - <Z 


(72) 

(78) 

(74) 

(75) 

(76) 

(77) 
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where 


£ Z = 


mm 

k,t 


o, 


ft 


Therefore from equations (68), (75), and (77), ||e^(i)|| is bounded by 

\\eF(t)\\<(T^-CWj(t) 


(78) 


(79) 


To bound ||e/f(i)|| apply equations (69), (75), and (79) to the norm of equation (66) and simplify to 
obtain 


)|ejc(0ll < (l + ^z) o^-CW 2 F {t) + ~h z oCW F (t ) 

Square equation (80) and use the fact that |WV(t)| < 1 to obtain 


(80) 


l|e/f(t)ll 2 < W 1 | (l + h) 2 ( J M-) Vj 

+ M 2 ( i+ *)(rHH (,) } 


(!) 


(81) 


The relative RMS error is defined as 


_ [/-^IIEkMH 2 ^] 172 [nje K (t)|| 2 dt 


e RMS = 


[roc l|F[<5W].||2d w ] 1/2 [/^ ||<5 W(f) || 2 dt 


1/2 


1/2 


(82) 


The inequality in equation (82) follows from Parseval's theorem (ref. 6) and inequality (67). The integral 
in the denominator of equation (82) is evaluated by means of equation (75) to obtain 


/ oo r oo 

||«5W(t)|| 2 dw = C 2 / Wj{t)dt 

- oo J — oo 


(83) 


Integration of equation (81) over t and division of the result by equation (83) yields the bound to 


'RMS 


as 


,2 <r „2 

e RMS - ° 


1 + 


’ll 2 

K 


)(r) 




1 + K 


where 


f-ooW^dt 

° = IIPaII • II Aj| 

hz = ||M ? || • ||P A II (j~^) \\F{d' 


^=-||Mj JP A || JF(6w)|| max 


2 0 + ^)(r) 4 ] 


e m — max 
k,t 


£z = 


min 

k.t 


o, 


6$k(ty 

ft 

ft 


(69) 

(58) 

(59) 
(73) 
(78) 
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and from equations (74) and (75) 


W P (t) = 11 v ;|1 

F[ j max||«W(f)|| 

The terms involving h z /K in equation (84) can be made negligible by scaling F(8 w) through suitable 
choice of factor K. Scaling F(6w) in effect limits the distortion caused by its projection from DFP 
coordinates to OBJ coordinates via mapping D. On the other hand, the error due to distortion caused 
by projection of 8W(t) from OBJ coordinates to DFP coordinates, represented by the remaining term in 
equation (84), is independent of scale factor K . This portion of the error is roughly proportional to the 
<; component of <5W(0, represented by the ratio 6^/(1 - c z ). 

If K is chosen sufficiently large so that the terms involving h z /K are negligible in equation (84), then 
f RMS i s bounded by 


f RMS < ° 



fZoWpWdt 


1/2 


1/2 


< a 


\ 

i -tz) 


(85) 


The second bound in equation (85) follows from the fact that since 0 < Wp{t) < I, then 


0 < wj(t) < wj(t) 


(86) 


and therefore the ratio of integrals in equation (85) is bounded by unity. 

Recall that ||A|| equals the largest singular value of A. It is easily shown that ||Pa|| equals the square 

root of the largest eigenvalue of ^A r A^j , which equals the inverse of the smallest nonzero singular 
value of A. Thus a equals the ratio of the largest singular value of A to the smallest nonzero singular 
value of A, defined as the condition number of A (ref. 5). Matrix A tends to be poorly conditioned (a 
large) if cameras 1 and 2 are closely spaced with nearly parallel focal planes (angles <j)\ and <f >2 small in 
fig. 4) or if they are nearly opposed (0 [ and (j ) 2 close to tt/2 in fig. 4). Matrix A is well conditioned (a 
close to 1) if the principal axes are coplanar and at right angles to each other ( (j)\ and <fi 2 equal to 7t/4 in 
fig. 4). 


Evaluation of Bound for Exponentially Damped Sinusoidal Perturbation Functions 

Inequality (85) is evaluated analytically for the special case where (*>W(/) is a vector of fixed direction 
of varying magnitude. Let 

FW{t)=hf{t) (87) 

where h is a vector in OBJ coordinates. 

A typical perturbation function encountered in structural vibration testing is the exponentially damped 
sinusoid 

/(<>={'"'“ m 

Tho minimum value of c~ ,lt sinu>o< occurs for u;q t approximately equal to 37r/2. Therefore 

(z = (M e- Sna / 2 ^ (89) 


Analytic evaluation of inequality (85) for this case yields the bound 


( RMS < a 


( ( M \ 

^0 ] 

U -<z) 

1 

(NO 

3 

+ 

CN 

00 


1/2 


( 90 ) 


A similar computation can be performed for the exponentially damped cosine perturbation function 


/(/) = 


(' COS u)[)t 
0 


(t > 0) 1 

(t < 0) / 


(91) 
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to give 


e RMS < ° 


{ \ ( 16u^ ■+■ 28 q^uJq -f- 3 u)q \ 

\ 1 — (-z ) y64o^ + 48 o 2 Wq + 8cjg J 


(92) 


For metal structures, typical values of a/uiQ seldom exceed 0.02, in which case inequalities (90) and (92) 
reduce to 


e RMS < 0.613 <t 



(93) 


Numerical studies showing computed errors for damped sinusoidal perturbations compared with maximum 
errors predicted by equation (93) are presented in reference 2. 


Geometric Examples 

Figure 4 illustrates an example of typical geometry in which the y\- and ?/ 2 -axes are parallel to the 
Z c -axis. and the coordinates (fi.fr) and (£ 2 ^ 2 ) are rotated from the (X c , Y c ) coordinates by angles 4>\ 
and 02, respectively. Let the object W be located at the intersection of the normals to each focal plane. 
Rotation matrices Mj and M 2 are then 


Mi. = 


cos sin 0*. 0 

0 0 1 

sin 0*. — cos 0*. 0 


(94) 


for k = 1,2. Matrices A and defined in equations (20) and (13), are constructed from Mi and M 2 

as 

'ci cos^i ci sin0i 0 ' 

0 0 ci 

C2 cos 02 ” c 2 sin 02 0 

- 0 0 c 2 . 


A = 


(95) 


and 


Mf = 


sin 0i —cos 0i 0 
sin 02 — cos 02 0 


(96) 


If constants are chosen such that 02 = — 0i, Cl = C2t ci = C 2 , 0i — 02 = tt/2, and 9\ = 0 2 = 0, then the 
matrix product A r A is 

cos 2 000 

0 sin 2 0 0 | (97) 

0 0 1 


A T A — 2 c 2 


The norm || A]| , equal to the square root of the largest eigenvalue of A T A (ref. 5), is then 

ii A II = y/2 c 


(98) 


The norm ||P A ||, equal to the square root of the largest eigenvalue of (A T A) 1 , is 



||P A || = max(sec 0, esc 0) 

(99) 

Similarly, it is determined that 

||M C || — \/2max(sin0,cos0) 

(100) 

For this case 

a — |[P A || ■ || A || = max(sec 0, esc 0) 

(101) 

Note that if <5W(<) is confined to the Z c -direction then 7i = 72 = 

7r/2, and hence by equation (73) 


e M = 0. 
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Camera 1 Camera 2 



Figure 5. Parallel projection geometry showing projection of object perturbation h W onto parallel focal planes. 


Iii the second geomet ric example, shown in figure 5, the OBJ and DFP coordinate systems are parallel; 
the camera focal planes are eoplanar; focal lengths rj and c 2 are equal; and points W Vv , W P2 , and W lie 
on the plane Y c — 0. Since <Hi/d <'<iuals Aq^/Ci for this case, it follows that 


AZ(Z + tfZ) _1 = 


? + 6$ 


Substitute equation (102) into equation (63), note that matrices and A cancel, and take the norm to 
obtain 




c + 0? 


aw) + ~h z \\F{6W)W 


Equation (85) simplifies to 


(■&■) [/“ wUt)dt] l/2 

CRMS < J ( 104 ) 

{fZcWpWdt} 

(o.quivalcnt to a — 1). Matrices M c and P^, necessary for determination of h z and v via equations (58) 
and (59) in the computation of scale factor K , are now obtained. Rotation matrix M is the identity 
matrix I, since the OBJ and FP coordinates are parallel. Matrix M ? is then 


and |)M s -j| = \/2. Matrix A is given by 


0 0 1 
0 0 1 


• 1 0 b 1 • 
0 1 0 
1 o b 2 
-0 1 0 . 
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where 


b Xi-Xp^ 


(*' = 1 , 2 ) 


Matrix A T A is 


A 1 A = 2 c 2 


1 o ^2 1 
0 1 0 
1^+6^ q 62+b| 


Its inverse is 


(A r A) -1 = 


(h - b 2 ) 2 c 2 

The eigenvalues of ( A r A)~ J are obtained analytically as 


62+62 

2 


0 


n (6i —62)^ 

h +&2 0 


0 

1 


A -J_ 1 

2c 2 ’ 2 


( a A + b\) ± J (a x + b\) 2 - 4 - df'j 


where 


bj + b 2 

°A = 7 - ,2 9 

(61 - 6 2 ) c 2 


b a = 


^A=- 


(61 - 6 2 ) 2 c 2 
6 l + 62 


(61 - 6 2 ) 2 c 2 

The norm of ||Pa|| is the square root of the largest eigenvalue of (A^A) -1 . If b\ — —b 2 , then 


H Pa| I = \/lc max(1, | ^ l|) 


(107) 

(108) 

(109) 

( 110 ) 

( 111 ) 

(112) 

(113) 

(114) 


Effects of Error on Estimated Frequency Response Functions 

The effects on estimated frequency response functions caused by the computation of F[<5W(£)] by the 
approximation KD F[<5w(f)] j can be deduced by examination of the error expression in equation (42), 
whose dominant part has been shown (see eqs. (63) and (70)) to be 

E /f ( W )=P A F[D,(f)AW(l)] (115) 

for suitably large K. Consider the jth element of D z (f) A<5W(f), which is 

< 116 > 

where A j. is the jth row of A and k = 1 when j — 1 or 2 and k = 2 when j — 3 or 4. Recall from 
equation (15) that 

«ft(t) = Mj(fiW(t) (117) 


19 


Expand (%/&) / [1 + (tfft/fc ) \ into a P owcr series 




which is convergent for |<5g c (£)/0c| < 1- For simplicity, suppose that SW(t) contains only two damped 
sinusoidal components at frequencies uq and u>2- Then 6${t)/<; will be of the form 

— = aie~ ait cosuqt + Oi^er a<li cosu^J (119) 


The jth component of A6W(t) will be of the form 

Aj. iW(() — f3\e~ (Llt cosuqZ + /? 2 C _a2 * cosu^ (120) 

From equations (116), (118), (119), and (120), a power series expansion for ej(t) is obtained as 

oo n 

ej{t) = Y2 b nm(t)cos n ~ m uiitc °s TO w 2 1 (121) 

n =2 to = 0 

where the coefficients b nrn (t) are weighted sums of products of powers of e~ a lt and e~ (l2t . By means of 
trigonometric identities one obtains 


m/2 

cos m (jjt — ^ c mq cos quit 
9=0 

where coefficients c rnq sat isfy 0 <'mq = T or 


(m even) 


(m — 1)/2 

cos m u )t = ^ c' mr/ cos quit 

Q = 0 

where coefficients satisfy c f mq ~ 1 and c m0 = 0, and finally 


(m odd) 


cos <juqJ cos = - [cos (guq + rcJ2) t — cos (guq - ru^) t] 
Substitute formulas (122) through (124) into equation (121), and simplify to obtain 


1 oc nr (n— m)/2 m/2 

«j(t) = oE E r nm ^ XI XI d n -m,qdm,T [<’°S (<7^1 + rw 2 ) i 
n—2 m=0 l g=0 r=0 


+ cos (qui - ruj 2 ) t] 


where d lJ = if i is even and dq = cF if i is odd. Equat ion (125) shows that the error frequency spectrum 
contains peaks at harmonic frequencies (all integer multiples of uq and u> 2 ) and at cross-modulated 
frequencies (sums and differences of all fundamental and harmonic frequencies). The amplitudes of the 
spectral peaks decrease rapidly with increasing frequency because the nth harmonic is weighted by (l/^) n 
in t he power series. 
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If the experimenter can identify true modes in D{F[w(t)]} and discard harmonic and cross-modulated 
error peaks, then error effects are relatively unimportant to modal analysis, since the fundamental 
frequencies are not shifted. However, estimated damping factors may be somewhat in error because 
of distortion in the shapes of fundamental peaks. 


Estimation of Frequency Response Matrix 

The dynamics of a three-dimensional structure, whose frequency response functions are to be identified, 
are modeled by a three-dimensional coupled system of iVth-order linear time-invariant differential 
equations in OBJ coordinates. Let 8XJ(t) be a three-element-vector driving function and let 8W(t) 
be a three-element-vector response function at point <5W. Then 6W{t) and 6U(t) are related by 

N d n 

£ ^ p n6W(t) + P 0 6W(t) = 6U(t) (126) 

71 = 1 


where Pq, Pi, 
obtain 


P n are 3x3 matrices of constants. Take the Fourier transform of equation (126) to 


N 


n = 0 


«W(w) =tfU(w) 


(127) 


By means of techniques of linear systems analysis (ref. 7), equation (127) may be solved for 8W(uj) as 


fiW(w) = H(w)«UM (128) 

where H(cj), defined as the 3x3 frequency response matrix, is the inverse of o(^) n Pn* H is desired 
to estimate H(cj) from observations in dual focal plane (DFP) coordinates. 

Let (eq. (73)) be sufficiently small, so that mappings D ^ and D , defined in equations (21) and 
(28), are approximated well by 


tfw(i) - D a [6W ( 0 ] « Z -1 A 6W{t) (129) 

and 

8W(t) — D[6w(t)] as PaZ 8w(t) (130) 

Although P^ is only a left inverse of A, it follows from equation (30) that 

Z" 1 APaZ <5w(t) = 6w(t) (131) 

whenever 8w(t) is contained in the column space of Z -1 A, or equivalently in the range of mapping D& 
Since all Sw(t) and 6u(t) to be considered are contained in the range of D a PaZ and Z* 1 A are hereafter 
treated as two-sided inverses of each other. Therefore 6JJ(t) and <5u(£) are related by 

6U(0«PaZ5u(0 (1.32) 


6u{t) w Z _1 AfiU(*) (133) 

Let 8JJ(t) be a zero-mean multivariate stochastic process which is stationary in the wide sense (ref. 8). 
Then <5W(£), 6u(i), and 8w(t) are also zero-mean wide-sense stationary processes. The cross covariance 
matrix between 6W(t) and 6XJ(t) is defined as 


Rwv{r)=E\ i 6W(t)8V* T {t-T) 


(134) 
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and the autoeovariance matrix of 6JJ(t) is defined as 


Ruu(r) = E [<5U(0<5U* T (t-r)' 

(135) 

where E denotes expected value. 


The cross-spectral matrix Swu(^) and the autospeetral matrix S 

lyu( w ) are defined as 

SwuM = ^PwuW] 

(136) 

SuuM = F [Ruu( r )l 

(137) 

Similarly define in DFP coordinates 


SuuM = F[E \bu(t) 6 u* T ( t - t) 

} (138) 

Swu(w) = [^w(0^u* T (< - t) 

} (139) 

It can be shown from equations (130) and (132) that 


SuuM = PaZSuuM zp^ 

(140) 

Swu( w ) = PaZ S wu (w) zp t a 

(141) 


It is shown in reference 9 that if />W(a>) and MJ(w) are related by equation (128), then Swuf^) and 
Syu(^) are related by 

SwuM — H(w) Syu(u;) (142) 


Thus, if S\vu( w ) an ( l SulK^) are known, the solution of equation (142) determines H(ce). It can 
be shown (ref. 9) that if W(£) contains additive noise, then the solution of equation (142) furnishes a 
miniinuni-inean-square-error estimate of H(u;). 

Since experimental data are recorded in DFP coordinates, it is advantageous computationally (ref. 2) 
to estimate S uu (w) and S W u(^) and to then estimate H(u;) by using equations (140) through (142). 
Smoothed est imates of S uu (oj) and S wu (^) are computed as averaged conjugate outer products of the 
Fourier transforms of N observed records of length T of Sv/(t) and 6u(t). Thus the Fourier transforms of 
6w(t) and fiu(t) for the nth record are 

^w n (u;) = [ 6w(t) e~ lujt dt (143) 

J (n— 1)T 

rnT 

6n{u)= 6u{t) (~ 1 ^ dt, 

J (n—\)T 

and the smoothed estimated autospeetral and cross-spectral matrices are 

Suu(w) = — (u) 

71 = 1 

1 N 

Swu(w) = — ^2 < >) w n (u;) (5u* r (cj) (146) 

n= 1 


(144) 


( 14 > r >) 
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Smoothed estimates of Suu(cj) and Swij(cj) may be obtained from S uu (o;) and S wu (u;) by means of 
transformations (140) and (141) to yield 

Suu( w ) = PaZ S U uM ZP a (147) 

SwuM = PaZS wu (w) ZP^ (148) 

In principle H(w) can be estimated by equation (142), provided Suu(w) is nonsingular, as 

H(w) = S wu (w) Sutj(w) (149) 

If the three components of U (£) are mutually independent stochastic processes, then Suu( a; ) is a 
nonsingular diagonal matrix, from which H (uj) is readily obtained via equation (149). If Su u(co) is 
nonsingular but nondiagonal, calculation of its inverse in equation (149) for each value of u may be 
excessively expensive computationally. If U (t) is confined to a line or to a plane in space, so that 

U(t) = Si(t)Ui+j2(t)U 2 (150) 

where Ui and U 2 are fixed, then Suu(cj) is singular, and equation (149) cannot be employed. Two 
alternative procedures for estimation of H(u;) are now proposed. 


Estimation of Diagonal Frequency Response Matrix 

The observed driving function, <5u(£), and the observed response function, 6w(t), in DFP coordinates 
are related through equations (126), (130), and (132). Take Fourier transforms of 6W(£) and 6TJ(t) in 
equations (130) and (132) and premultiply equation (128) by Z _1 A to obtain 


<5w(cj) = Z~ 1 AH(o;)PaZ(5u((j) (151) 

It follows from equation (151) that S W u(^) and S uu (u;) are related by 

Swu(w) = Z _1 A H(u-’) PaZS uu (w) (152) 

so that. 

PaZ Swu( U i1 ) — H(co’) PaZ S uu (cj) (153) 

If it is known that H(u;) is diagonal (interactions between the X c , Y c , and Z c components are absent), 
then H (oj) is easily obtained from equation (153) as follows. Let 

^wu(^) — PaZ S W U M (154) 

and 

Suu(^) = PaZSuuM (155) 

so that equation (153) becomes 

SwuM = H(cj) S uu (cj) (156) 

Then h^(uj) can be estimated by using equation (156) as 


^tu( w ) 


^uu^M 


(157) 


for i = 1, 2 , 3, and j = 1, 2, 3, 4, where <r wu tj (uj) and a uu ^ ( cj) are elements of X) wu (u;) and £ uu (cj) 

and * denotes values obtained from the estimated matrices S wu (^) and S uu (w). If H( cj) is diagonal, 
then hnj(u) will be independent of j. Conversely, if equation (157) is independent of j for arbitrary 
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S wu (^) and Euu(w), then H(w) must be diagonal. Thus, the diagonality of H(w) may be confirmed by 
computations of h U j(^) (eq. (157)) for all 12 combinations of indices for i = 1, 2. 3 and j = 1, 2, 3, 4 
and verification that, h^ j(oj) is independent of j. 


Estimation of Nondiagonal Frequency Response Matrix 

In the general case where H(u;) is a full nondiagonal matrix, let a unidirectional driving function g(t) 
bo applied in three independent measurements in the X c -. Y c -, and ^-directions separately. One column 
of H(w) may then be estimated from each measurement., in DFP coordinates, of the observed input <5 u ( / ) 
and the observed response <5w(f). Thus apply g(t) in the Ac-direct ion. Input vector 5U(£) is then 


SV{t) = 9(t) 


1 

0 

0 


Expressed in DFP coordinates, the input vector is 


6u(t) = g(t)Z~ l A 


1 

0 

0 


Autospectral matrix S uu (cj) is then 


Suu(w) 


,s '(/!/(w)Z A 


1 0 
0 0 
0 0 


0 

0 

0 


A r Z“‘ 


(158) 


(159) 


( 160 ) 


where .s (/!/ (w) is the autopower spectrum of g(l). E<]iiation (160) may be rewritten as 

P A Z S UU (W) = ,Syy(cj) Tfc 

where (for k = 1) 


r, = 


"11 "21 "31 "41 

Cl Cl C2 C'2 

0 0 0 0 

0 0 0 0 


( 161 ) 


(162) 


and djj is the clement in the ?‘th row and jth column of A. From equations (153) and (101), it follows 
that 


p A zs wu (U) = Sgy(w)H(u)r k 


(103) 


Antopower spectrum s f/r/ (cj) can be estimated from S U u(^) by means of equation (101). Column 1 of H(cj) 
is then estimated from the estimated s^(cj) and cross-spectral matrix S wu (u;) by means of equation (163). 

Columns 2 and 3 of H(u>) are estimated in a similar manner. Apply g(t) in the V^-direction for A; = 2 
to give 


and 


*U (/) = g{t) 


0 

1 

0 


r 2 = 


'±12 

Cl 

0 


0 

"22 

Cl 

0 


0 0 1 


"32 

C2 

0 


"42 

C2 

0 


> 


(164) 


(165) 
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Use equations (161) and (163) to estimate s gg { cj) and column 2 of H(u>). 
Z c -direction for k = 3 to obtain 


*U(t)=ff(t) 


0 

0 

1 


Finally apply g(t) in the 


(166) 


and 


r 3 


0 0 0 O' 

0 0 0 0 

a 13 023 0.33 043 

?1 ?1 ?2 ?2 - 


(167) 


Use equations (161) and (163) as before to estimate s gg (oj) and column 3 of H(c j). 


Least-Squares Estimation of Frequency Response Matrix 

A method is developed for least-squares estimation of H(cj) from equations (161) and (163). Since 
Suu(^) and S wu (u;) are estimated by means of equations (145) and (146) and since measurements of 
u (t) and w(£) are noisy, no values of s ffff (u;) and H(cj) exist in general which will satisfy equations (161) 
and (163) for estimated matrices S uu (oj) and S wu (u;). Thus, in general for every value of s gg (u) in 
equation (161) 

Suu(^) 7 ^ Sggi 1 ^) (168) 

where 

A/s = Z _1 Ar fc (169) 

and for every H(cj) in equation (163) 

P A zS wu ( W )^S w ( W )H( W )r fc (170) 

Optimum estimates of s gg (uj) and column k of H(cj), denoted by h.^(u;), can be determined by choosing 
s gg {uj) and h./^u;) such that the sum of the squares of the Euclidian distances between corresponding 
columns is minimized in inequalities (168) and (170). Thus, from equation (168) define 

4 

e u ~ ^2 “ Sggi^) l| 2 

3~ 1 
4 

= ^2 [®UU.j ( W ) _ Sgy[(d) \.j 

3 = 1 


§uu.j (w) - 


(171) 


where s uu .(w) is the j th column of S uu (u), and A.y is the j th column of A*. 

To minimize cj, differentiate equation (171) with respect to s gg {u), set the result equal to zero, and 
solve for s gg (u) to obtain 


>'. 4 


To estimate h./^u;), use equation (170) to define 


d = E II^U.,-(w) - 8gg{u) Ik] h.j(w)|| 2 

3 = 1 


(173) 
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where is the j th column of PaZS wu (^) and s gg (u) is given by equation (172). Proceed as 

above to minimize and solve for h.^(w) to obtain 




lkj° WU.J (^) 

*w(") 


where 7^ is the element in the fcth row and jth column of matrix r*.. 


(174) 


Concluding Remarks 

The accuracy of frequency response function estimation from stereoscopic dynamic data by a rapid 
computational procedure has been computed. It was shown that the error incurred by reversal of 
the operations of coordinate transformation and Fourier transformation includes terms dependent on 
the magnitudes of components of motion normal to the camera focal planes (whose effects cannot be 
eliminated), as well as on the condition of the stereoscopic projection matrix. Additional error terms exist 
which may be eliminated by proper scaling of the Fourier transforms prior to coordinate transformation. 
Spectral analysis of the dominant error term, whic h is dependent on the scale factor, shows the presence 
of harmonic and cross-modulated spectral peaks caused by the nonlinearity of the perspective projection 
from spatial coordinates to focal plane coordinates. Only when the camera geometry is such that all 
motion is parallel to both camera planes is the spatial transformation linear and no error introduced. 

The full 3x3 frequency response matrix of a three-dimensional structure can be practicably estimated 
from one stereoscopic dynamic data set only if the orthogonal components of the driving function are 
mutually independent, stochastic processes. Alternative procedures have been given for estimation of a 
diagonal frequency response matrix and for least -squares estimation of a full matrix, colunm-by-column 
in three independent sets of measurements. 
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